Prom In-In to Pb-Pb collisions: Lessons from high-precision Dilepton Measurements 
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The NA60 experiment has studied low mass muon pair production in 158 AGeV In-In collisions 
with unprecedented precision. With these measurements one can constrain the in-medium modifica- 
tions of hadrons more thoroughly than from previous experiments. We present a model description 
which is able to describe the data well. We argue that model constraints could be significantly 
further tightened and theoretical uncertainties in modeling of the medium evolution reduced by an 
additional experiment measuring low mass di-muon production in Pb-Pb collisions at 158 AGeV 
with the same resolution as provided by NA60, thus allowing to disentagle the in-medium spectral 
function from effects of the evolution model with greater precision. 

PACS numbers: 25.75.-q 



I. INTRODUCTION 

The changes of properties of vector mesons (most promi- 
nently the p^) in a hot and dense medium have been 
studied extensively in model calculations since their elec- 
tromagnetic decay into lepton pairs makes them experi- 
mentally accessible. Hadronic in-medium properties have 
also been linked to the chiral phase transition, since vec- 
tor and axial- vector correlators are predicted to merge as 
the chirally restored phase is approached. However, the 
ideas about how this is realized in detail diverge. Brown 
and Rho proposed that the restoration of chiral 

symmetry might manifest itself in a scaling of the vec- 
tor meson masses with the order parameter of the phase 
transition. On the other hand hadronic many-body cal- 
culations differ in the predicted mass shifts but agree on 
the conclusion that the p-meson spectral function is sig- 
nificantly broadened due to its interactions with the sur- 
rounding medium [1 H 0, II H 113 ■ Early QCD sum rule 
calculations seemed to validate dropping mass scenarios 
[H 113, but QCD sum rules have later been shown to 
constrain the correlation between in-medium width and 
mass and not on one of the quantities alone 0, 0| . Ul- 
timately, detailed knowledge of all hadronic sources of 
dileptons is also crucial to address the question if in high 
energy nucleus-nucleus collisions new non-hadronic con- 
tributions to the measured dilepton spectrum arise, such 
as radiation from dcconfincd matter. 
The CERES and HELlOS-3 collaborations have shown 
that central nucleus- nucleus collisions exhibit a strong 
enhancement of low mass dileptons compared to scaled 
proton-nucleus and proton-proton collisions |l5l IT^ 
Hq . Although these measurements demonstrate that in- 
deed considerable in-medium modifications of the vector 
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meson channel occur, their precision is not high enough 
to conclusively discriminate between a wide range of dif- 
ferent scenarios 0, 0, ^ . Not even the rather straight- 
forward question whether a mass shift of the p-meson 
due to the Brown-Rho scaling law Q, |^ H, Q is realized 
in nature or not can be answered conclusively Q from 
these results. 

The high precision di-muon data of the NA60 experiment 
now provide the possibility to put strong constraints 
on our theoretical understanding of in-medium modifi- 
cations of vector mesons. Most model calculations for 
in-medium modifications so far have made the simplify- 
ing assumption that the p-meson is at rest. Thus, the 
only reasonable comparison between calculations using 
this assumptions and the data employs the NA60 results 
at low transverse lepton pair momentum. However, since 
this data set is not yet published we postpone a system- 
atic and extensive theoretical study comparing different 
in- medium models to a forthcoming publication [l9l | . We 
will nevertheless report some of the results of this study 
[T^ here. We find that the NA60 data require consider- 
able in-medium broadening of the p-meson while strong 
mass shifts arc excluded. Moderate mass shifts are not 
favored by the NA60 data and can be clearly ruled out if 
not accompanied by strong broadening. 

In the present paper we follow a different idea and ar- 
gue that the spectral function favoured by the data in 
combination with an evolution model based on a scaled 
version of a well-constrained Pb-Pb evolution model is 
able to provide a good qualitative and quantitative de- 
scription of the NA60 data in all momentum bins. We 
then turn with the same spectral function to the Pb-Pb 
system where the evolution model is much better con- 
strained by hadronic data po| and show what additional 
insights beyond what is seen in In-In could be gained 
if the Pb-Pb system could be measured with the same 
precision as In-In. 
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II. CALCULATING THE DIMUON SPECTRUM 

Schematically, the emitted spectrum of dileptons from 
an evolving thermalized system can be found from the 
convolution 



dMdpTdr] 



evolution ( 
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—r— (E) acceptance (1) 
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where the rate is found from the spectral function 
R{q,T,pB) as 



dN 



a" R{q,T,pB) 



(2) 



Here, q is the 4-momentum of the emitted muon pair, 
T is the temperature of the emitting volume element, 
13 = 1/T and pB its baryon density. The fireball evolu- 
tion encodes information on radiating volume, T, baryon 
chemical potential pB, transverse flow velocity vt, longi- 
tudinal rapidity rj and chemical non-equilibrium proper- 
ties such as pion chemical potential /i^. 
In the case of NA60, the acceptance is highly non-trivial 
(and for the comparison with the measured spectrum, a 
full Monte Carlo simulation of the detector has been car- 
ried out by NA60 with the calculation as input). Thus, 
while a Lorentz-invariant quantity like the invariant mass 
spectrum dN/dM is in principle independent from the 
flow profile of the medium, the resulting invariant mass 
spectrum after acceptance folding is not: The acceptance 
has pt dependence and thus the hardening of the momen- 
tum spectra due to flow alters the shape of the invariant 
mass distribution. Thus, for a meaningful comparison of 
theoretical spectral functions with the data it is manda- 
tory that both the medium evolution and the acceptance 
are implemented in more than just a schematic way. 



A. The evolution model 

A detailed description of the model used to describe the 
fireball evolution is found in po| . The main assumption 
for the model is that an equilibrated system is formed a 
short time tq after the onset of the collision. Further- 
more, we assume that this thermal fireball subsequently 
expands isentropically until the mean free path of par- 
ticles exceeds (at a timescale r/) the dimensions of the 
system and particles move without significant interaction 
to the detector. 

For the entropy density at a given proper time we make 
the ansatz 



siT,7^„r)=NR{r,T)-H{7j„T) 



(3) 



with T the proper time as measured in a frame co-moving 
with a given volume element, r the radius, 77^ the space- 
time rapidity and i?(r, r), Hijjs, t) two functions describ- 
ing the shape of the distribution and a normalization 



factor. We use Woods-Saxon distributions with a r de- 
pendent scale Rc{t) {Hc(t)) and a thickness parameter 
dws iVws) to describe the shapes for a given r. For a 
radially non-relativistic expansion and constant acceler- 
ation we find Rc{t) ~ Rq + ■ Hc{t) is obtained 
by integrating forward in r a trajectory originating from 
the collision center which is characterized by a rapidity 
rjcir) = ?7o + a-TiT with rjc = atanh where is the 
longitudinal velocity of that trajectory. Since the rela- 
tion between proper time as measured in the co-moving 
frame and lab time is determined by the rapidity at a 
given time, the resulting integral is in general non-trivial 
and solved numerically (see [23| for details). Rq is deter- 
mined in overlap calculations using Glauber theory, the 
initial size of the rapidity interval occupied by the fire- 
ball matter. 770 is a free parameter and we choose to use 
the transverse velocity = o,±Tf and rapidity at decou- 
pling proper time rj-^ = ijq + a^^Tf as parameters. Thus, 
specifying 770, ?7/, and r/ sets the essential scales of the 
spacetime evolution and dws and 77ws specify the detailed 
distribution of entropy density. 

For transverse flow we assume a linear relation between 
radius r and transverse rapidity p = atanh v±{t) = 
r/Rc{T) ■ Pc{t) with Pc{t) = atanh a±T. We allow for the 
possibility of accelerated longitudinal expansion which in 
general implies rj ^ rjs pol |. 

The model parameters have been adjusted to hadronic 
transverse mass spectra and HBT correlation measure- 
ments in 158 AGeV Pb-Pb and Pb-Au colhsions at SPS. 
The framework was then used to successfully describe 
photon and dilcpton emission and charmonium suppres- 
sion |20l |. Since at this point no freeze-out analysis or 
HBT correlation data are available from NA60 for the 
In-In system, our strategy is to start with the parame- 
ter set derived for Pb-Pb collisions and use geometrical 
scaling arguments to go to In-In. 

We infer the change in total entropy production from 
the ratio of charged particle rapidity densities dNch/drj 
at 30% peripheral Pb-Au collisions with 2.1 < ?/ < 
2.55 measured by CERES [H El El and semi-central 
In-In collisions at 77 = 3.8 measured by NA60 [2ll |. 



dN: 



-In 



Pb-Pb 



0.68. The number of participant 



baryons and the effective initial radius (for the sake of 
simplicity, we map the overlap area in non-central colli- 
sions to a circle with the same area) are obtained with 
nuclear overlap calculations. We assume that the stop- 
ping power (which determines the width of the inital 
distribution of entropy and baryon number in rapidity) 
scales approximately with the number of binary collision 
per participant, Nhin/Npart- This determines rjQ. The 
shape parameters d^s and rj^js have no great impact on 
electromagnetic emission into the midrapidity slice, they 
primarily govern the ratio of surface to volume emission 
of hadrons. Thus, we leave them unchanged from their 
value determined in |20l | . Under the assumption that the 
physics leading to equilibration is primarily a function of 
incident energy we keep the formation time r = 1 fm/c 
as in Pb-Pb collisions. We stress that the final results 
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exhibit no great sensitivity to the choice of cither rjo or 
To except in the high M, high px hmit. 

The biggest uncertaintly is the choice of the decouphng 
temperature Tp. Due to the fact that the In-In system 
is smaller than Pb-Pb, we expect a higher decoupling 
temperature. However, without simultaneous knowledge 
of HBT correlations and transverse mass spectra an un- 
ambiguous answer cannot be obtained. We tentatively 
choose Tp = 130 MeV. Our lack of detailed knowledge 
translates primarily into an uncertainty in the normal- 
ization of the final results (which depend on the emission 
duration which increases for lower Tp). It is for this 
reason that the result is normalized to the data in the 
interval M < 0.9 GeV in the comparison to account for 
discrepancies of the order of a few percent. 
We fix the flow velocity at Tc by choosing the value of 
transverse flow obtained at T = 130 MeV in our Pb-Pb 
model in the same centrality class. There is likewise an 
uncertainty (which is to some degree correlated with the 
choice of Tp) since flow determines the slope of hadronic 
(and to a lesser degree electromagnetic) px spectra. 
The equation of state in the hadronic phase and the off- 
equilibrium parameters /i7r,/iK are inferred from statis- 
tical model calculations as described in . The re- 
sulting fireball evolution is characterized by a peak tem- 
perature of about 250 MeV, a lifetime of ^ 7.5 fm/c and 
a top transverse flow velocity of 0.5c at decoupling. 



B. The spectral function 

Below the phase transition temperature the effective de- 
grees of freedom are hadrons. The photon couples now 
to the lowest-lying dipole excitations of the vacuum, the 
p, iu and (j) mesons and multi-pionic states with the ap- 
propriate quantum number. 

There arc different effective hadronic models and tech- 
niques for the calculation of the properties of hadronic 
matter near and below the phase boundary. For reviews 
see e.g. |^ |3- Most of the effective model calcu- 
lations that are constrained by phenomenological data 
such as the hadronic and electromagnetic decay widths 
of the particles predict substantial broadening in matter 
with comparably small shifts of the in-medium masses. 
As already discussed in the introduction we focus here 
on one specific hadronic model calculation which is 
favoured by the data and discuss it as a 'generic' exam- 
ple. The model calculation employed here solves trun- 
cated Schwindcr-Dyson equations incorporating a self- 
consistent resummation of the tt— p interaction in order to 
determine the thermal broadening of the p-meson. For a 
detailed discussion of this approach see [lOj and citations 
therein 12^. This ^-derivable approach self-consistently 
takes into account the finite in-medium damping width 
of the pion in the thermal heat bath which contributes 
substantially to the broadening of the p-mcson. Effects of 
baryons are not included in this approximation scheme. 



We note that the finding that only thermal effects would 
be able to account for the data is clearly interesting. 
As long as the thermodynamically active degrees of free- 
dom are quarks and gluons, the timelike photon couples 
to the continuum of thermally excited qq states and sub- 
sequently converts into a charged lepton pair. The cal- 
culation of the photon spectral function at the one-loop 
level is performed using standard thermal field theory 
methods. The leading-order result for bare quarks and 
gluons as degrees of freedom is well known j^- This re- 
sult, however, holds only up to perturbative higher order 
corrections in gs that take into account collective plasma 
effects. We describe those by means of a quasiparticle 
model, additional details can be found in 0. 

C. The vacuum p- meson 

In addition to the dilepton radiation from the medium 
(mainly through the p'^ channel) there is also radiation 
from e.m. decays of p mesons in vacuum after the ther- 
malized evolution has ceased. In principle there is also 
radiation from co and (f> vaccum decays, however due to 
the small width of these states this contribution is exper- 
imentally very well identified and has been subtracted 
|2l| . The vacuum decay of p mesons is commonly re- 
ferred to as 'cocktail p' and obtained by calculating the 
produced number of p mesons by a statistical model cal- 
culation. 

In such a calculation, the number of p mesons is found 
under the assumption that the relative abundances of 
hadron species are given by the thermal expectation value 
above a temperature Tc (close to the transition temper- 
ature Tp) but that subsequently only resonance decays 
occur. In such a framework, the decay of heavy reso- 
nances leads to an overpopulation of pion phase space 
as compared to the pion phase space density in chemical 
equilibrium. This can be parametrized by the introduc- 
tion of a pion chemical potential p^ (which grows with 
decreasing temperature) and leads via the tightly coupled 
p — > TTTrp channel to an effective fip = 2^-^. 
However, while tracking of the resonance decays fixes the 
number of p-mesons at thermal decoupling, their momen- 
tum spectrum remains unknown. In particular, p mesons 
after chemical decoupling also participate in the collec- 
tive expansion of matter. While this effect does not in- 
fluence the integrated invariant mass spectrum, the ac- 
ceptance leads to a significant influence of flow on the de- 
tected signal. Thus, we go beyond the statistical model 
formulation and compute the full momentum spectrum of 
p mesons after decoupling using the Cooper- Frye formula 

^^ = (2^y^^^^^^^piHr^J 

with the momentum of the emitted p and da^ an 
element of the freeze-out hypersurface (determined by 
the condition T = Tp) just as we calculate the emission of 
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any other hadron from the medium. We let this ensemble 
undergo electromagnetic decay with vacuum width and 
lifetime. 



III. RESULTS 

In Fig. ^ left panel, we present the comparison of the 
calculation to the NA60 all pT data for In-In collisions 
at 158 AGeV using the spectral function obtained in the 
<f>- functional approach [lOj . The calculation contains the 
vacuum pas described in section III CI and a DD back- 
ground |2l|. We note that the agreement with the data 
indicates that the relative contribution of the vacuum 
p and in-medium p is at least approximately verified. 
While the strength of the vacuum p parametrically scales 
like the freeze-out volume times the p lifetime V{Tf)TP^^, 
the in-medium radiation scales like J V{T)dT and is thus 
sensitive to the lifetime of the fireball. Thus, the relative 
strength is a consistency test for the medium lifetime in 
the In-In collision evolution model. 

Fig. n E^lso presents the calculation of the di-muon spec- 
trum that would be expected in 158 AGeV Pb-Pb colli- 
sions with the same detector and resolution as used by 
the NA60 collaboration in In-In collisions. In order to 
compare In-In collisions and Pb-Pb collisions we normal- 
ized both calculation to the In-In collision data in the 
mass region below 0.9 GeV. We note that the model de- 
scribes the absolute normalization of the In-In data set 
with less than 5% discrepancy. 

In Fig. ^ right panel, we compare the predictions for the 
Pb-Pb system and In-In system for low pT normalized 
to NA60 In-In low pr-data (the yet unpublished low pr 
data are not shown) . This shows that the low pT 
cut acceptance of the detector puts serious limitations 
to the capability of low px data to distinguish different 
lineshapes, in particular in the mass range below the p 
mass. 

Comparing both calculations one realizes the following 
most prominent distinctions between the spectra for Pb- 
Pb collisions and In-In collisions: 

1. The relative contribution of the in- medium p to the 
vacuum p is more pronounced due to the longer 
lifetime in Pb-Pb. 

2. Contributions at low invariant mass are more em- 
phasized in the Pb-Pb system due to lower Tp. 
This is clearly seen in the all pr binning but some- 
what obscured for the low pr data. 

These features are a clear test for different evolution mod- 
els and in particular for differences in lifetime and Tp. 
They show that the in-medium modifications of the p- 
meson would show up more prominently in a 158 AGeV 
Pb-Pb measurement with the same resolution than in 
In-In, since the relative contribution of the vacuum p- 
meson would be reduced. The Pb-Pb collision experi- 
ment would therefore probe the in-medium modifications 
in the mass region around and below the vacuum rho- 
mass more clearly than in the In-In system at all and 



high pt- However, it would also probe these spectral 
functions at on average lower temperature. This could 
provide information of how the the in-medium properties 
vary in the region between the different Tp of the Pb-Pb 
and In-In systems, going beyond what can be constrained 
by In-In alone. 

A third distinction, though expected, is not visible in the 
results. Since the radiation from hot volume elements 
with T > Tc is increasingly important for the contri- 
bution to invariant masses above the vacuum p-mass, 
it is therefore expected that the higher initial temper- 
ature of the Pb-Pb system leads to a harder slope of the 
Pb-Pb di-muon spectrum at high invariant masses where 
more relative contribution to the total yield is expected 
from above T > Tc than in In-In. Although this is true 
when comparing the spectra before acceptance, after ac- 
ceptance correction the effect is not significant in the all 
Pt data. 

IV. CONCLUSIONS AND OUTLOOK 

In this paper we studied the di-muon production in 
158 AGeV In-In collisions measured by the NA60 collab- 
oration. We first presented a model which reproduces the 
In-In data qualitatively and quantitatively well. We then 
argued that a measurement of Pb-Pb with the same accu- 
racy would be capable of providing additional non-trivial 
information and so help to disentangle effects of spectral 
function and evolution model. Already the pT-mtegrated 
lineshape indicated sizeable differences which in the sim- 
ulation arise from different lifetimes and freeze-out tem- 
peratures. It is of considerable advantage that the dy- 
namics of Pb-Pb collisions is already well constrained 
from measurements of liadronic observables. The expec- 
tation in the high mass region would also differ in com- 
parison between Pb-Pb and In-In collisions. We argued 
that this could provide important information about con- 
tributions from non-hadronic sources. This is especially 
important since a purely hadronic evolution seems un- 
able to reproduce the NA60 In-In data in the high pT 
high M region which is an unexpected and interest- 
ing hint that even in the rather small In-In system some 
amount of partonic evolution may contribute to the ob- 
served radiation. This contribution should be increased 
in Pb-Pb collisions and be visible in the Pb-Pb to In-In 
comparison if a high pT binning is performed. 
We will present a detailed systematic and extensive theo- 
retical study comparing different in-medium models and 
their predictions for different pt cuts for the In-In system 
in the near future Il9l. 
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M < 0.9 GeV. 
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